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Energy distribution in nonaxisymmetric magnetic Taylor- Couette flow 
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Azimuthal magnetorotational instability is a mechanism that generates nonaxisymmetric field pattern. Nonlinear simula- 
tions in an infinite Taylor-Couette system with current-free external field show, that not only the linearly unstable mode 
m = 1 appears, but also an inverse cascade transporting energy into the axisymmetric field is possible. By varying the 
Reynolds number of the flow and the Hartmann number for the magnetic field, we find that the ratio between axisymmetric 
(m = 0) and dominating nonaxisymmetric mode (m = 1) can be nearly free chosen. On the surface of the outer cylinder 
this mode distribution appears similarly, but with weaker axisymmetric fields. We do not find significant differences in the 
case that a constant current within the flow is added. 
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1 Introduction 

MHD-instabilities in radiation zones of differentially ro- 
tating stars influence the rotation profile of the star due 
to angular momentum transport (Denissenkov & Pinson- 
neault 2007; Yoon, Langer & Norman 2006). Tayler insta- 
bility (Tayler 1973) or helical magnetorotational instability 
(Hollerbach & Riidiger 2005) are examples for such insta- 
bilities. According to Spruit (1999) the Tayler instability 
could possibly support a dynamo. In their model the in- 
stability should create a significant poloidal field that due 
to differential rotation is wound up into a toroidal one and 
thus close the loop for a self-sustained dynamo (Braithwaite 
& Spruit 2004; Braithwaite 2006). While this idea sounds 
promising, there are doubts about its general functionality. 
Simulations in very similar geometry with z-dependent dif- 
ferential rotation did not show a dynamo (Gellert, Riidiger 
& Elstner 2007). Zahn, Brun & Mathis (2007) claim that 
such a process is not possible in this simple way. The au- 
thors argue that the nonaxisymmetric modes wound up by 
differential rotation can not contribute to the axisymmetric 
field, because the azimuthal wavenumber is always m 7^ 0. 
As we show in the following, azimuthal magnetorotational 
instability (AMRI) creates also an axisymmetric component 
in the nonlinear regime. For a scenario based on AMRI as 
instability mechanism that might overcome the arguments 
of the incompatible wavenumber and would make field re- 
generation at least possible. It does not necessarily mean, 
however, that it supports a dynamo. 

This work consists of two aspects. First we describe the 
development of nonaxisymmetric modes in the nonlinear 



regime of AMRI-unstable Taylor-Couette flow. Especially 
the observed inverse cascade transporting energy into the 
axisymmetric field and its dependence on the shear is an- 
alyzed. As second aspect we point out that weak currents 
within the flow do not basically influence the instability. 
Field structure and energy distribution remain almost un- 
changed. 

2 Equations and numerical treatment 

We use the hydrodynamic Fourier spectral element code de- 
scribed by Fournier et al. (2005) extended by its magnetic 
field. With this approach we solve the 3D MHD equations 
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(rotB) X B, (1) 



dtu + (it • V)u = -Vp + V^u 

dtB = —V'^B + rot(M X B), (2) 
Pm 

divM = 0, divB = (3) 
for an incompressible medium in cylindrical coordinates 
(i?, 0, z). Free parameters are the Hartmann number 
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Ha = S 

and the magnetic Prandtl number Pm = v/rj. Here v is 
the viscosity of the fluid and ?/ its magnetic diffusivity. The 
Reynolds number is defined as Re = VliaRi^D / v with D ~ 
Rout — Rin (unit of length) and the angular velocity of the 
inner cylinder 57in. Unit of velocity is ly/D and unit of time 
the viscous time /v. 

The solution is expanded in M Fourier modes in the 
azimuthal direction. This gives rise to a collection of merid- 
ional problems, each of which is solved using a Legendre 
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spectral element method (see e.g. Deville, Fischer & Mund 
2002). Either M = 8 or Af = 16 Fourier modes are used, 
three elements in radial and eighteen elements in axial di- 
rection. The polynomial order is varied between = 8 
and N = 12. With a semi-implicit approach consisting of 
second-order backward differentiation formula and third or- 
der Adams-Bashforth for the nonlinear forcing terms time 
stepping is done with second-order accuracy (Fournier et al. 
2004). 

At the inner and outer wall perfect conducting boundary 
conditions are applied. In axial direction we use periodic 
boundary conditions to avoid perturbing effects from soUd 
end caps. The periodicity in z is set to F = QD. 

The initial flow profile is the typical Couette profile 

n{R) =a+^ (5) 
with 

a = r^i^in, b = r^i? Siin, (6) 

V = -Rout/^in = 0.5 and fiQ = fiout/^^in = 0.5. This 
means the outer cylinder rotates with half of the angular ve- 
locity of the inner cylinder in the same direction. After the 
Rayleigh stability criterion dR{R^ft)^ > this configura- 
tion is hydrodynamically stable. In addition a toroidal exter- 
nal field -Bext — (0, -B°, 0) is applied. The second compo- 
nent has a radial profile similar to the flow; 

B°iR)=aBR+^. (7) 
The external field than can be characterized by the number 

Bout asRout + bsRout .OS, 

MB - -Ho- - n • 

-Din flB^in + Ofiilin 

The current-free field is given by fiB = 0.5. Values be- 
low lead to a field connected with constant negative cur- 
rent jz within the flow and larger values to a positive cur- 
rent. To trigger the instability the initial magnetic field con- 
sists of small random perturbations (with an amplitude of 

3 Results 

3.1 Linear instability 

From linear stability theory is known that only the mode 
TO = 1 is unstable (Rudiger et al. 2007a), which is reflected 
in the nonlinear simulations, too. Only m = 1 starts to grow 
in the beginning of a simulation where one is in a quasi- 
linear regime. The stability boundaries found with the non- 
Unear code compared to linear results are shown in Fig.[T]for 
Pm = 1,/ta = 0.5 and fis — 0.5. Equal to the results in 
Rudiger et al. (2007b) we find not only an onset of the insta- 
bility if the shear (given by the Reynolds number) exceeds 
a certain threshold for a given magnetic field configuration 
(constant Ha), but also an upper boundary. Strong shear pre- 
vents the instability - from this point of view an instability 



AMRI (instability of m=1) for Pm=1 , n,j=0.5 




60 80 100 120 140 160 180 200 
Ha 

Fig. 1 Stability diagram for Pm = l,fin = 0.5 and 
fiB — 0.5 (AMRI, current-free external field). The solid line 
results from linear stability analysis, dots and stars mark re- 
sults from nonlinear simulations. The dashed line is an as- 
sumed continuation of the solid line, there were no linear 
data available. 

of the magnetic field. Both, the linear results and nonlinear 
simulations, agree very well on the lower boundary. At the 
upper boundary for Re > 300 the simulations exhibit a shift 
of the stability boundary to slightly higher values of Re. But 
the qualitative picture is the same. 

Following a line with constant Re, one can argue that 
weak but not to small magnetic fields lead to an instability 
of the flow. For higher Re also the magnetic field has to 
exceed a higher threshold to sustain the nonaxisymmetric 
mode generation. Insofar both descriptions, unstable flow 
due to a magnetic field and magnetic field instability due to 
shear, are equivalent. 

3.2 Nonlinear simulations 
3.2.1 Current-free external field 

After the onset of the instabihty, first exclusively the to = 1 
mode grows consistent with the linear theory. If the to = 1 
mode reaches roughly 30% the amplitude of the final state, 
also both neighboring modes to — 0, 2 start to grow. 
Slightly later successively higher and higher modes follow. 
The energy drops by a factor of ten for each to m + 1. 
Modes TO > 3 carry less than lO^'^ of the total energy and 
are not plotted in Fig. |4] A time series of the total magnetic 
energy 

E{m)^^ j^B\m)AV (9) 

for each of the significant modes is shown in Fig.|2l the flow 
and magnetic field pattern is shown in Fig. |3] Compared 
to the strength of the external field, the magnetic field due 
to the AMRI is weaker. The maximum size is around 0.3 
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Fig. 2 Time series of the total magnetic energy for Fourier 
modes < m < 3 for Rc = 200. Linearly unstable is only 
m = 1. Time is scaled in periods of rotation. 
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Fig. 3 On the left an isosurface plot of 
(±0.7|BK|max = 0.09BP„) for Re = 200, Ha = 110 
at t = 300 is shown, on the right the magnetic field 
component in the i?-z -plane also in units of BjJ^. 



BP„ and the energy below 10%. The saturated state shows 
no regular time-dependence, the field pattern rotates with a 
small drift relative to the system rotation (as already pointed 
out by the linear analysis). 

Interesting is the fact that in the nonlinear simulations 
the modes m > 1 appear only very weak. The energy of 
m = 5 is akeady by a factor of 10^ less than the energy 
of the mode to = 1. On the other hand also a significant 
large mode to = develops. This means energy transport 
happens not only to higher modes or smaller scales but also 
in the opposite direction 0. 

By varying the shear strength, the amount of energy 
in the produced axisymmetric mode can be nearly free 
chosen. Near the stability boundary (e.g. Re = 175 and 



' Regarding a dynamo based on the instability, Zahn et al. (2006) point 
out that the generation of an axisymmetric field from a nonaxisymmefiic 
one is difficult. If the instability itself also creates an axisymmetric mode 
the generation of an appropriate poloidal axisymmetric field might be eas- 
ier to achieve. 






Fig. 4 Total magnetic energy distribution in Fourier 
modes to with growing shear for Ha ~ 110. The instability 
is observable between Re = 175 and Re — 390. 



Re = 390 for Ha — 110, see Fig. [TJ the nonaxisymmet- 
ric mode m — 1 dominates the energy spectrum, in be- 
tween at Re = 280 to = becomes of the same size. 
This seems to be a general behavior and is not only the 
case for fin — 0.5. Also for steeper rotation laws, e.g. the 
quasi-keplerian (/tQ = 0.3535), this kind of energy trans- 
port into the axisymmetric mode is observable in the simu- 
lations. The amount of energy in the axisymmetric mode is 
always of the order of 50% of the whole energy or less. This 
agrees with magnetic fields observed in Ap stars. Accord- 
ing to this observations slow rotating Ap stars with strong 
fields often possess a nonaxisymmetic field. And the axi- 
and nonaxisymmetric modes are of the same size (Land- 
street & Mathys 2000; Bychkov, Bychkova & Madej 2005). 
Regarding these observations more interesting than the vol- 
ume averaged energy is the behavior of the magnetic field 
that penetrates the surface and would be visible from out- 
side. This is the poloidal part of the field, which consists 
due to the boundary conditions only of the z-component. 
Looking on the poloidal field, the situation is similar to the 
global field distribution. The dominating mode is yet the 
TO = 1 mode. Different is the behavior of the produced ax- 
isymmetric mode. It is nearly not visible on the surface. And 
opposite to the global behavior it becomes even less with 
Reynolds numbers in the middle of the unstable parameter 
space, where the energy of the to = mode is maximal in 
the volume-averaged quantity. This is illustrated in Fig. |5] 
where the power spectral density \Bz{m)\'^ at R — i?out 
and z = l/3r is plotted for different Reynolds numbers. 
The created axisymmetric component contributes mainly to 
the i?0-component, the original source of the instability. 
Due to the perfect conducting boundary conditions, which 
are quite different to the vacuum environment for real stars, 
the surface field is not directly comparable with the obser- 
vations. 
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Fig. 5 Power spectral density for on the cylinder sur- 
face for several m normalized by the total power on the sur- 
face. The distribution is similar to the volume-averaged val- 
ues, except the m = mode occurs to be smaller in Bz- 
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Fig. 6 Influence of a positive or negative current on the 
energy distribution (jiB — 0.35,0.5,0.65 from left to 
right) is not visible. The axisymmetric mode caiTies Eo = 
0.249, 0.265, 0.272 % of the total energy, the m = 1 mode 
El = 0.690, 0.685, 0.682 % for /is = 0.35, 0.5, 0.65. His- 
tograms are normahzed with total magnetic energy. 



4 Conclusion 

In an infinite Taylor-Couette system we have shown, that 
nonlinear simulations of AMRI lead to an unstable nonax- 
isymmetric mode m = 1 of the magnetic field if a certain 
strength of the external magnetic field is reached. Further- 
more the instability is suppressed if the external field be- 
comes too strong (for given Reynolds number). Equivalent 
is the point of view of an instability of the magnetic field 
due to shear. Here also a certain strength of the shear or 
value of the Reynolds number has to be reached to trigger 
the instability of the to = 1 mode. And if the shear exceeds 
a certain value, the instability is not observable anymore. 

The simulations reveal in the saturated state a distribu- 
tion of magnetic energy over several modes, not only the 
ni — 1 is excited. However, the energy transport works 
in both directions, not only in the small-scale direction 
to larger m, but also to lower to. A significant amount 
of energy is shifted inversely into the axisymmetric mode 
(m = 0). The shear strength is responsible for the ratio 
between magnetic energy in axisymmetric and nonaxisym- 
metric mode TO = 1. Near the stability boundaries, the non- 
axisymmetric mode caiTies nearly all energy. In the central 
region of the unstable parameter space the energy ratio be- 
tween modes to = and to = 1 is roughly one. 

Adding a small positive or negative current within the 
flow does not affect the instability, the field pattern and final 
energy distribution do not change compared to the current- 
free case. A small deviation from the profile cx 
does not visibly influence the occurrence of the instability. 
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3.2.2 External field connected with a current 

With the profile of the external field as defined in Eq. (|7]i, 
an associated current is always constant. It is positive if 
fiB > 0.5 and negative if p,B < 0.5. To test the influence 
of its current on the appearance of the instability and the 
resulting energy distribution, we alter the external field pro- 
file X.O fiB — 0.35 and /is = 0.65. These modifications 
of the external field are small enough not to change the pro- 
file qualitatively. The simulations with the modified external 
fields do not show a basic influence in the energy distribu- 
tion from a cuiTent. The ratio between magnetic energy of 
modes to = and m = 1 has values of Eq/Ei = 0.37 
for fiB = 0.35, £^0/^1 0.39 for fiB = 0.50 and 
Eq/Ei = 0.40 for fiB ~ 0.65. With increasing jiB also 
the energy of the external field increases and is responsi- 
ble for the slightly increasing energy ratio. The instability is 
not strongly effected by low currents. For larger currents, 
or larger Hartmann numbers, the influence of the Tayler 
instability grows and the Reynolds number dependence of 
the onset of the instability is lost. Tayler instability already 
works with low or even without shear. 
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